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EXECUTIVE  SUMMARY 


The  purpose  of  this  project  was  to  refine  and  further  develop  an  instrumented 
mannequin/flash  fire  exposure  system  developed  in  an  earlier  project,  and  to  use  the  system  to 
study  the  effect  of  several  variables  on  the  measured  protective  qualities  of  various  protective 
garments.  This  report  describes  the  results  obtained. 

The  number  of  skin  simulant  sensors  in  the  mannequin  was  increased  from  80  to  110,  and 
the  computer  controlled  data  acquisition  system  modified  to  accept  the  additional  30  input 
channels.  In  addition,  the  computer  codes  were  modified  to  record  and  store  the  additional  data. 
To  produce  more  accurate  skin  bum  predictions,  the  thermal  conductivity  of  the  skin  simulant 
material  was  measured  with  a special  instrument  at  the  National  Research  Council  of  Canada. 
The  measured  value  was  about  100%  larger  than  the  value  used  previously. 

The  mannequin/flash  fire  system  was  installed  in  a new  facility  dedicated  for  this  purpose. 
Additional  burners  and  gas  storage  capacity  were  added  to  the  flash  fire  generating  system.  This 
permitted  steadier  gas  flows  and  more  uniform  heat  flux  exposure  to  the  mannequin  surface.  The 
repeatability  in  the  predicted  surface  area  of  the  mannequin  receiving  second  degree  bums  or 
worse  is  better  than  95%  for  average  heat  fluxes  up  to  84  kW/m2  and  up  to  4 seconds  duration. 

Testing  of  representative  garments  showed  that  the  effects  of  varying  heat  flux  and 
exposure  times  differs  with  garment  type,  and  that  while  the  TPP  test  results  showed  similar 
trends,  the  TPP  results  did  not  predict  the  behaviour  of  all  of  the  garments  tested  on  the 
mannequin.  It  is  suggested  that  whenever  possible,  the  exposure  conditions  for  mannequin 
testing  should  be  selected  to  most  closely  simulate  the  hazardous  environment  form  which 
protection  is  sought. 
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A.  INTRODUCTION 


In  many  industrial  settings  workers  potentially  can  be  exposed  to  flash  fires  and  fireballs. 
The  oil  and  gas  and  petrochemical  industries  are  ones  where  ignition  of  gas  leaks  can  produce 
accidents  resulting  in  injury  and  loss  of  life  through  skin  bums.  Investigation  of  numerous 
accidents  suggest  that  workers  can  be  exposed  to  intense  heat  flux  (=  80  kW/m2)  for  relatively 
short  periods  of  time,  typically  five  seconds  or  less.  While  protective  clothing  is  available  for  use 
in  these  potentially  hazardous  environments,  evaluation  of  such  clothing  typically  has  been  through 
small  scale  laboratory  tests  of  flame  resistance  and  thermal  protective  performance.  Although  such 
tests  may  give  valuable  information  about  the  fabrics  used  in  thermal  protective  clothing,  they  do 
not  necessarily  predict  how  well  garments  or  garment  assemblies  will  perform  when  actually 
exposed  to  a flash  fire. 

One  method  of  evaluating  protective  garments  is  to  expose  a clothed  instrumented 
mannequin  to  a controlled  flash  fire  and  to  assess  the  resulting  potential  skin  damage.  The 
potential  skin  damage  can  be  evaluated  using  computer  models  if  the  local  heat  flux  to  the  surface 
is  known.  If  the  heat  flux  and  exposure  time  can  be  varied  reliably,  then  different  accident 
scenarios  can  be  simulated.  The  problem  addressed  by  this  research  was  the  refinement  and 
further  development  of  an  instrumented  mannequin/flash  fire  exposure  system  developed  in  an 
earlier  project,  "Systematic  Assessment  of  Protective  Clothing  for  Alberta  Workers"  [1].  This  is 
the  only  mannequin  of  its  type  in  Canada  and  one  of  the  few  world-wide.  The  purpose  of  the 
system  is  to  assist  industry  in  assessing  the  protective  qualities  of  garment  systems  when  subjected 
to  flash  fires  of  varying  duration  and  intensity,  by  estimating  potential  for  skin  bum  damage. 

The  system  as  originally  developed  was  limited  in  its  uniformity  of  exposure  and  range  of 
heat  flux.  In  addition,  the  estimation  of  potential  for  skin  damage  is  directly  dependent  upon 
knowing  the  thermal-physical  properties  of  the  skin  simulant  sensors.  Thus,  accurate  measurement 
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of  these  properties  is  required  to  obtain  the  best  possible  estimates.  The  purpose  of  this  research 
was  (a)  to  develop  and  refine  an  instrumented  mannequin  and  flash  fire  exposure  system,  and  (b) 
to  determine  the  effect  of  varying  heat  flux  and  exposure  time  on  the  measured  protective 
properties  of  several  protective  garments.  Specific  objectives  were: 

(a)  to  modify  the  flash  fire  generation  system  to  provide  both  (i)  controlled  variation  of  heat 
flux,  and  (ii)  uniform  exposure; 

(b)  to  measure  the  thermal-physical  properties  of  the  skin  simulant  material; 

(c)  to  add  more  skin  simulant  sensors  to  more  thoroughly  cover  the  mannequin  as  well  as 
additional  data  acquisition  capability; 

(d)  to  modify  the  computer  code  to  speed  calculations  and  provide  better  graphical  portrayal 
of  skin  damaged  areas;  and 

(e)  to  test  various  garment  systems  at  varied  heat  fluxes  and  exposure  times. 

The  remainder  of  the  report  is  divided  into  two  main  sections.  The  first  section  covers 
objectives  (a)  to  (d),  relating  to  the  refinement  of  the  mannequin  and  flash  fire  exposure  systems. 
The  second  section  reports  on  the  experiments  in  which  heat  flux  and  exposure  time  were  varied. 
A brief  section  on  interlaboratory  testing  follows.  A final  section  includes  a summary  of  the 
project  and  our  conclusions. 

B.  EXPERIMENTAL  SYSTEM  AND  METHOD  OF  EVALUATION 

The  male  mannequin,  size  40R,  was  constructed  out  of  fibreglass  by  duplicating  the  shape 
of  an  existing  store  mannequin.  This  was  done  using  female-male  molding  techniques.  Hand  lay 
up  of  matting  and  resin  was  used  to  produce  a 3 mm  thick  shell.  This  thickness  of  fibreglass  was 
used  in  order  to  provide  sufficient  strength  for  handling  and  to  have  a thick  enough  layer  for  easy 
mounting  of  heat  flux  sensors.  Before  construction,  samples  of  hardened  fibreglass  were  exposed 
to  large  flash  fires  at  a fire  fighting  school  to  ensure  suitability. 
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Figure  1 shows  the  front  and  rear  views  of  the  mannequin.  A lift  out  panel  was  designed  into 
the  chest  area  to  facilitate  the  installation  of  the  heat  flux  sensors  and  their  wiring.  The  output 
cable  from  all  the  sensors  is  fed  out  of  the  mannequin  through  a ten  cm  diameter  aluminum  pipe 
molded  into  the  head.  A steel  cable  wrapped  around  this  pipe  is  used  to  hang  the  mannequin. 
Note  that  the  mannequin  exhibits  large  muscle  groups  typical  of  those  found  on  a well  conditioned 
adult  male.  The  arms  and  hands  are  removable  to  facilitate  dressing.  Mounting  sockets  were 
molded  into  the  shoulders  and  forearms  for  this  purpose.  The  hands  are  solid  fibreglass. 


Fig.  1 Front  and  rear  views  of  instrumented  mannequin 
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Heat  Flux  Sensors 

The  thermal  protective  quality  of  a garment  is  judged  on  estimates  of  the  extent  of  skin  damage 
resulting  from  a controlled  flash  fire.  In  order  to  estimate  the  skin  damage,  the  heat  transfer  to 
the  surface  of  the  mannequin  needs  to  be  measured.  Once  the  heat  transfer  and  its  variation  with 
location  and  time  is  known,  a heat  transfer  model  of  human  skin  can  be  used  to  estimate  the 
damage  due  to  burning.  The  heat  flux  sensors  (skin  simulant  sensors)  used  in  the  mannequin 
should  ideally  be  as  close  as  possible  to  human  skin  in  their  response  to  the  heat  transfer.  This 
minimizes  the  corrections  that  have  to  be  applied  in  the  model. 

Heat  transfer  from  a short  duration  flash  fire  into  the  surface  of  the  mannequin  (or  the  skin  of 
a human)  can  be  modelled  as  transient  heat  transfer  into  a semi-infinite  slab.  For  the  case  of  a 
constant  heat  flux  at  the  surface  and  constant  thermal  physical  properties  the  analytical  solution 
can  be  shown  to  be  [2] 


T(x,t)  -T:  = 1 

' k 

Where: 

T = temperature,  K 
Tj  = initial  temperature  in  slab,  K 
q = heat  flux,  W/m2 
k = thermal  conductivity,  W/m-  K 
a = thermal  diffusivity  (k/pc),  m2/s 
p = density,  kg/m3 
c = specific  heat,  J/kg-°C 
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e = 2.7183 

x = distance  into  slab,  m 
erfc  = complimentary  error  function 
At  the  surface  (i.e.  x = 0)  this  reduces  to 

r(o,f)  - Ti  = (2) 

yjnkpc 

This  latter  form  illustrates  that  the  product  kpc,  the  thermal  inertia,  or  its  square  root,  the  thermal 
absorptivity,  is  the  grouping  of  thermal  physical  properties  that  should  be  matched  between  the 
heat  flux  sensors  and  human  skin. 

The  design  selected  for  the  skin  simulant  sensors  is  shown  in  Fig.  2.  The  material  "Colorceran" 
is  commonly  used  to  form  the  top  layer  of  chemistry  laboratory  benches.  The  thermal  physical 
properties  of  this  material,  at  least  as  far  as  its  ability  to  simulate  the  short  duration  heat  transfer 
from  a flash  fire  into  human  skin  is  concerned,  are  close  to  ideal.  Table  1 lists  values  of  the 
thermal  physical  properties  of  the  two  outer  layers  of  human  skin,  that  is  the  epidermis  and  dermis 
layers,  and  those  of  Colorceran.  While  the  individual  properties  are  not  a perfect  match  with  those 
of  the  epidermis  and  dermis  layers,  the  product  kpc,  or  its  square  root,  are  closely  matched.  The 
thermal  physical  properties  for  human  skin  listed  in  Table  1 were  taken  from  the  literature  [3,4]. 

The  shape  of  the  skin  simulant  sensors  was  selected  for  convenience.  The  Colorceran  comes 
in  32  mm  thick  slabs.  The  19  mm  diameter  was  selected  to  give  a large  enough  diameter  so  that 
sufficient  lead  length  could  be  given  to  the  thermocouple  on  the  surface  to  reduce  conduction  heat 
transfer  effects  to  the  junction.  A hole  was  drilled  along  the  long  axis  of  the  sensor  to  allow  the 
thermocouple  wire  to  be  run  into  the  inside  of  the  mannequin.  The  thermocouple  wire  used  was 
30  gauge  copper  constantan.  The  section  on  the  surface  of  the  sensor  was  rolled  flat. 
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approximately  0=  1 mm  thick,  prior  to  installation.  The  wire  is  held  on  to  the  surface  of  the  sensor 
with  a special  high  temperature  epoxy-phenolic  adhesive.  The  adhesive  has  a maximum  short  term 
working  temperature  of  370°C,  which  was  thought  to  be  well  above  the  surface  temperatures 
reached  in  the  short  duration  flash  fires.  To  date  no  sensor  failures  have  occurred  due  to  the 
adhesive  losing  its  bond. 


30  gauge  Cu-Co  wire,  rolled  flat 
~ 0.1  mm  thick  held  on  surface 
with  an  epoxy-phenolic  adhesive. 

"Colerceran"  - inorganic  mixture  of 
calcium,  aluminum,  silicate  with 
asbestos  fibres  and  a binder. 


Fig.  2 Skin  simulant  sensor 
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TABLE  1 Thermal  Properties  of  Human  Skin  and  Skin  Simulant 


Property 

Human  Skin 

Skin 

Simulant 

Epidermis 

Dermis 

k (W/m-  K) 

0.255 

0.523 

0.97 

p (kg/m3) 

1200 

1200 

1877 

c (J/kg-K) 

3598 

3222 

1205 

kpc  (J2/m4  • “C2  • s) 

LI  x 106 

2.0  x 106 

2.2  x 106 

Vkpc  (J/m2  • °C  • s1/2) 

1050 

1414 

1483 

To  adequately  cover  the  mannequin,  110  sensors  were  used.  Each  sensor  was  mounted  as 
flush  as  possible  with  the  fibreglass  exterior  of  the  mannequin  by  drilling  a hole  in  the  shell  and 
gluing  the  sensor  in  place  with  fibreglass  resin.  The  sensors  can  be  removed  by  cutting  the  resin 
bond.  The  distribution  of  the  sensors  is  shown  in  Fig.  3,  approximately  one  per  unit  surface  area. 
None  are  mounted  in  the  lower  arms  (wrist  region),  hands  or  feet  below  the  ankles.  After  the 
sensors  were  installed,  the  entire  mannequin  surface  was  painted  with  flat  black  high  temperature 
paint.  This  paint  has  a rated  temperature  of  600°C. 

The  thermocouple  wire  from  each  sensor  is  positioned  inside  the  mannequin  to  terminate  on 
a juncton  board  located  in  the  chest  cavity.  About  15  m of  copper  wire  is  used  to  connect  each 
thermocouple  to  the  data  acquisition  system.  The  extension  wire  exits  the  mannequin  through  the 
aluminum  tube  molded  into  its  head.  An  aluminum  foil  radiation  shield  is  used  to  protect  the 
exposed  extension  wire  from  the  flash  fires. 

Each  sensor  was  calibrated  in  situ  using  a 380  W movie  projector  lamp  powered  by  a 
variable  voltage  power  supply.  The  lamp  output  was  calibrated  with  three  slug  calorimeters  so  that 
its  radiant  heat  flux  was  known  to  within  a few  percent.  However,  the  calibration  accuracy  with 
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Fig.  3 Sensor  areas  on  instrumented  mannequin 
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the  movie  projector  lamp  is  limited  because  it  does  not  produce  a perfectly  uniform  heat  flux  on 
the  surface  of  the  sensor. 

The  thermal  absorptivity  of  each  sensor  was  determined  from  data  collected  during  the 
calibration  process.  The  specific  value  measured  for  each  sensor  was  stored  in  the  computer 
controlled  data  acquisition  system  and  is  used  whenever  heat  transfer  data  is  collected  during  a 
flash  fire.  Calibration  of  the  110  sensors  showed  that  the  mean  value  of  the  thermal  absorptivity 
was  within  7%  of  the  value  listed  in  Table  1,  and  that  the  standard  deviation  for  the  sample  was 
less  than  10%  of  the  mean.  A histogram  of  the  sensor  calibration  is  included  in  the  Appendix. 
The  thermal  absorptivity  value  listed  in  Table  1 was  calculated  from  measured  thermal  physical 
properties  of  the  sensor  material.  These  measurements  were  made  by  the  Institute  for  Research 
in  Construction,  National  Research  Council  of  Canada. 

FLASH  FIRE  GENERATION  SYSTEM  AND  FIRING  ROOM 

The  flash  fires  are  generated  with  propane  diffusion  flames.  Six  identical  burner  assemblies, 
positioned  around  the  mannequin,  are  used.  Each  burner  assembly.  Fig.  4,  consists  of  two  "tiger 
torches"  supplied  with  propane  at  three  atmospheres  absolute  pressure.  The  air  mixing  nozzles  in 
the  burner  heads  were  removed  and  flow  control  orifices  installed  upstream.  Interchangeable 
orifices  are  used  to  set  the  mass  flow  rate  of  fuel  through  each  burner.  The  three  atmosphere  gas 
pressure  was  selected  to  produce  choked  flow  in  the  orifices.  This  permits  the  average  heat  flux 
to  be  varied  from  about  67  kW/m2  to  84  kW/m2  (1.6  cal/cm2  • sec  to  2.0  cal/cm2*  sec).  Electric 
solenoid  valves  are  used  to  turn  on  and  off  the  gas  flow  to  the  burners.  The  duration  of  the  bum 
is  programmed  through  the  data  acquisition  system.  Pilot  lights  are  u§ed  to  ignite  the  main  gas 
flow  exiting  each  burner. 
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Fig.  4 Flash  fire  burner  system 
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The  average  heat  flux  was  determined  by  direct  measurement.  The  nude  mannequin  was 
subjected  to  either  a three  or  four  second  flash  fire  and  for  each  skin  simulant  sensor  the 
instantaneous  heat  flux  calculated  at  one  second  intervals.  The  mean,  maximum,  minimum  and 
standard  deviation  were  then  calculated  from  the  data  set.  The  mean  was  then  taken  to  be  the 
average  heat  flux  when  testing  garments  at  the  same  settings.  To  illustrate  values,  a typical  3 
second  exposure  would  produce  maximum  and  minimum  values  of  160  kW/m2  and  48  kW/m2, 
respectively  with  a mean  of  79  kW/m  and  a standard  deviation  of  22  kW/m  based  on  all  110 
sensors. 

The  flash  firing  system  is  housed  in  a special  concrete  block  room  added  on  to  the 
Mechanical  Engineering  Acoustics  and  Noise  Unit,  a facility  operated  by  the  Department  of 
Mechanical  Engineering,  University  of  Alberta.  The  room,  Fig.  5,  is  approximately  4 m by  4 m 
in  plan,  and  4 m high.  It  is  fitted  with  two  remote  controlled  dampers  for  fresh  air  supply  during 
firing  and  venting.  The  size  of  each  damper  is  approximately  0.5  m x 0.75  m.  Venting  of  gases 
is  done  after  firing  via  a fan  located  in  the  ceiling  of  the  room.  A 1.5  m by  1 m viewing  window 
made  of  wire  reinforced  glass  is  located  in  the  wall  that  separates  the  firing  room  from  the  main 
laboratory  building.  As  indicated  in  Fig.  5,  the  six  burner  assemblies  are  located  in  the  room  so 
as  to  produce  as  uniform  heat  flux  on  the  surface  of  the  mannequin  as  possible.  The  six  burner 
assemblies  are  fed  from  four  0. 1 m gaseous  propane  cylinders  set  in  the  comers  of  the  room.  The 
four  tanks  are  connected  by  steel  pipe  to  one  another  and  to  a 0.1  m3  liquid  propane  storage  bottle 
located  outside  the  room.  This  latter  bottle  is  fitted  with  a pressure  regulator  for  setting  the  gas 
storage  pressure  in  the  firing  room. 

The  burner  system  is  prevented  from  accidental  firing  through  safety  lockouts.  These  include 
a thermocouple  based  pilot  flame  sensor  system,  one  for  each  burner  head,  and  a normally  closed 
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solenoid  valve  used  to  isolate  the  outside  liquid  propane  storage  system  from  the  inside  vapour 
storage  containers.  Only  when  the  outside  tank  is  isolated  and  all  the  pilot  lights  are  burning  will 
the  computer-data  acquisition  system  permit  a signal  to  be  sent  to  start  the  gas  flow  to  the  burner 
heads. 

With  the  system  as  shown  in  Fig.  5,  only  a 1 1%  drop  in  mass  flow  rate  occurs  during  a three 
second  firing.  A ten  second  firing  produces  a 25%  decrease  in  mass  flow  rate  from  beginning  to 
end. 


Fig.  5 Layout  of  flash  fire  room 
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DATA  ACQUISITION  SYSTEM 

A computer  controlled  data  acquisition  system  is  used  to  run  the  tests,  record  and  store  data, 
and  to  calculate  and  display  the  results  of  the  numerical  model  used  to  estimate  the  skin  damage. 
The  output  from  the  110  thermocouples  on  the  skin  simulant  sensors  is  first  fed  into  an  analog 
amplifier  multiplexer,  the  output  being  amplified  200  times.  This  amplified  output  is  then  fed  into 
a 12  bit  analogue-to-digital  convertor  mounted  inside  the  computer.  The  digitized  signals  are 
written  either  directly  into  a temporary  memory  or  permanent  storage  by  writing  directly  to  a disk. 
The  conversion  rate  for  the  analogue  to  digital  converter  is  up  to  800  Hz,  thus  each  of  the  110 
thermocouples  can  be  read  about  eight  times  a second.  With  a 12  bit  system  the  signal  resolution 
is  1 part  in  4096,  while  the  800  Hz  sampling  rate  produces  a timing  error  of  about  0.125  seconds 
for  any  thermocouple. 

To  run  an  experiment,  software  specially  written  for  this  purpose  is  first  loaded  into  the 
computer.  It  prompts  the  user  for  information  about  the  test  such  as  the  bum  duration  and  total 
data  sample  time.  Once  this  is  keyed  in,  the  computer  system  checks  that  the  pilot  lights  for  the 
burners  are  lit  and  waits  for  a start  signal.  After  receiving  a start  signal  the  unit  displays  all  110 
temperatures  of  the  sensors.  If  these  are  satisfactory,  for  example  not  more  than  25 °C,  a run  signal 
is  given  by  the  operator.  The  system  then  takes  zero  readings  from  all  1 10  sensors  and  turns  on 
the  burners  for  the  prescribed  time.  Sampling  during  firing  is  done  at  800  Hz  during  burner  on 
times  and  then  continues  at  110  Hz  for  a maximum  of  two  minutes.  Once  sampling  is  finished 
the  temperature  data  is  loaded  into  a file  for  later  use. 
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SKIN  DAMAGE  EVALUATION 

The  purpose  of  the  protective  clothing  is  to  prevent  skin  damage  when  exposed  to  a flash 
fire.  Skin  damage  is  usually  identified  as  being  first,  second  or  third  degree  burn.  For  this  study 
only  second  degree  or  worse  is  of  importance. 

Many  investigators  have  examined  skin  burning  and  developed  criteria  to  predict  when 
irreversible  damage  begins  [5-9].  Damage  begins  when  the  growing  or  basal  layer  has  its 
temperature  raised  above  44°C.  At  this  temperature  irreversible  damage  begins  to  occur  to  the 
basal  layer  and  the  longer  it  is  at  or  above  this  temperature  the  greater  the  amount  of  damage. 
Experiments  on  human  and  pig  skin  showed  that  the  rate  of  destruction  of  the  basal  layer  could 
be  adequately  described  as  a chemical  reaction  of  the  first  order.  The  criterion  used  was  suggested 
originally  by  Henriques  [6]  and  is  of  the  form 


_A E 

BL  = Pe  RT  (3) 

dt 

where 


£2  = a quantitative  measure  of  bum  damage  at  the  basal  layer  or  at  any  depth  in  the  dermis 
P = frequency  factor,  s'1 
AE  = the  activation  energy  for  skin,  J/mole 
R = the  universal  gas  constant,  8.315  J/kmol*K 

T = the  absolute  temperature  at  the  basal  layer  or  at  any  depth  in  the  dermis,  K 
t = total  time  for  which  T is  above  44°C  (317.15  K) 

The  total  burn  damage  is  found  by  integrating  equation  (3)  over  the  time  interval  that  the 


basal  layer  is  above  44°C.  That  is, 
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i _ A E 

Q = J Pe  RT  dt  (4) 

o 

Henriques  [6]  found  that  if  Q is  less  than  or  equal  to  0.5  no  damage  will  occur  at  the  basal 
layer.  If  Q is  between  0.5  and  1.0  first  degree  burns  will  occur,  while  if  Q > 1.0  second  degree 
burns  will  result  This  total  damage  criteria  can  be  applied  to  any  depth  of  skin  provided  the 
appropriate  values  of  P and  AE  are  used. 

Morse  et  al.  [7]  have  examined  the  various  values  of  P and  AE  available  in  the  literature  and 
suggests  that  the  criteria  developed  by  Stoll  and  co-workers  [8,9]  be  used  in  the  epidermal  layer 
and  that  of  Takata  [10]  be  used  in  the  dermal  and  subcutaneous  layers.  The  values  of  P and  AE 
developed  by  Stoll  are: 


P = 2.185  x 10124  sec’1  - 
A E/R  = 93,534.9  K 

P = 1.823  x 1051  sec’1  • 
A E/R  = 39,109.8  K 

while  those  of  Takata  are: 


for  T < 50  °C 


for  T > 50°C 


P = 4.32  x 1064  sec"1  ► 
A E/R  = 50,000  K 

P = 9.39  x 10104  sec"1  > 
A E/R  = 80,000  K 


for  T < 50°C 


for  T > 50 °C 


In  order  to  predict  the  extent  of  damage  that  would  result  in  a flash  fire  it  is  necessary  to 
know  the  temperature  history  of  the  basal  layer.  The  heat  flux  sensors  in  the  mannequin  are  used 
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for  this  purpose.  They  permit  the  surface  heat  flux  and  its  variation  with  time  to  be  predicted. 
This  information  is  then  used  in  a numerical  model  for  heat  transfer  in  human  skin  in  order  to 
predict  the  temperature  history  of  the  basal  layer  and  the  extent  of  skin  damage  for  each  of  the  110 
sensors.  A full  description  of  skin  bum  criteria  can  be  found  in  the  Appendix. 

HEAT  TRANSFER  MODEL  OF  HUMAN  SKIN 

The  skin  model  used  is  based  on  one  developed  by  Metha  and  Wong  [3].  This  model  divides 
the  skin  into  three  tissue  layers,  the  epidermis,  dermis,  and  subcutaneous,  plus  an  isothermal  core. 
It  also  allows  for  blood  flow  through  arteries  and  veins  to  the  epidermis-dermis  interface.  An 
energy  balance  on  a region  of  skin  containing  both  arteries,  veins  and  capillaries  produces, 

cp  = k^L  - ( cp)bG(T-Tb ) (5) 

M dx2 

where: 

(cp)b  = volumetric  heat  capacity  of  the  blood,  J/m3  • K 
G = blood  perfusion  rate,  m /s/m  of  tissue  volume 

Tb  = blood  perfusion  temperature,  K (assumed  equal  to  the  internal  body  temperature, 

310  K) 

T = tissue  temperature,  K 

and  the  other  terms  are  as  defined  earlier. 

This  equation  is  solved  numerically  for  each  of  the  1 10  sensors.  The  object  is  to  predict  the 
temperature  distribution  and  history  in  the  skin  until  no  further  damage  occurs,  or  the  entire 
thickness  of  the  epidermis  and  dermis  layers  has  an  £2  value  greater  than  1.  For  the  heat  fluxes 
and  durations  used  in  this  study  the  numerical  model  always  predicted  that  no  further  skin  damage 
was  occurring  after  one  minute  after  the  flash  fire. 
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The  heat  flux  at  the  surface  of  the  mannequin  and  its  variation  with  time  is  determined  from 
the  recorded  temperature  data  from  the  sensors.  Duhamel’s  Integral  Method  is  used  to  convert  the 
measured  surface  temperature  of  the  sensor  into  a surface  heat  flux.  The  integration  is  done 
numerically  by  the  trapezoidal  rule.  This  heat  flux  is  then  used  as  the  surface  condition  in  the  skin 
heat  transfer  model. 

The  numerical  model  employs  varying  thermal  physical  properties  with  depth  as  shown  in 
Table  2.  The  Crank-Nicolson  Implicit  method  is  used  with  a stretching  ratio  for  the  grid  spacing 
in  each  of  the  three  layers.  The  grid  spacings  selected  were  six  in  the  epidermis,  24  in  the  dermis 
and  15  in  the  subcutaneous.  Several  different  time  steps  and  stretching  factors  were  tested  to  see 
which  produced  the  closest  result  to  a base  case  which  could  be  calculated  by  hand.  The  final 
form  uses  one  second  time  steps  and  a stretching  ratio  of  1.25. 


TABLE  2 Parameters  Used  in  Numerical  Model  of  Human  Skin 


Skin  Location 

Parameter 

Epidermis 

Dermis 

Subcutaneous 

Tissue 

Blood 

Overall  thickness  of 
layer  (m) 

8 x 10'5 

2 x 10'3 

1 x 10'2 

... 

Thermal  conductivity, 
k (W/m  • K) 

0.255 

0.523 

0.167 

... 

Volumetric  heat 
capacity,  cp  (J/m3  • K) 

4.32  x 106 

3.87  x 106 

2.76  x 106 

3.99  x 106 

Blood  perfusion  rate,  G 
m3/s/nr  of  tissue 

... 

1.25  x 10'3 

1.25  x 10‘3 

As  a test  of  the  skin  simulant  sensor  - skin  burn  - heat  transfer  model,  the  calibration  lamp 
was  used  to  heat  one  of  the  sensors  on  the  mannequin.  The  resulting  data  was  loaded  into  the  heat 
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transfer  model  in  order  to  predict  the  time  required  for  a second  degree  bum  to  occur.  Five 
different  heat  fluxes  were  used.  These  ranged  from  5.9  kW/m2  to  31  kW/m2  (0.14  cal/cm2*  s to 
0.74  cal/cm2’  s).  The  results  are  shown  in  Fig.  6.  The  time  to  achieve  a second  degree  burn 
agrees  very  well  with  the  experimental  results  of  Stoll  and  Quanta  [9].  This  indicates  that  the 
mannequin  system  should  give  reasonable  predictions  of  the  extent  and  nature  of  skin  burns  in  the 
flash  fire  simulations. 

The  details  of  the  computer  code  and  the  sensor  and  skin  models  are  given  in  the  Appendix. 
Included  are  descriptions  of  how  the  numerical  code  was  tested  and  overall  calibration  of  the 
sensor-skin-bum  model. 


Fig.  6 


Human  skin  tolerance  time  to  absorbed  energy  showing  calibration  points  obtained 
using  a calibrated  heat  flux  source  on  a skin  simulant  sensor 
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C.  GARMENT  TESTING  WITH  VARIED  EXPOSURES 
Six  different  protective  garment  types  were  exposed  to  simulated  flash  fires  varying  in  heat 
flux  and  exposure  time.  Exposures  of  67,  75  and  84  kW/m2  (1.6,  1.8  and  2.0  cal/cm2*  s)  for  three 
seconds  and  67  and  75  kW/m2  for  four  seconds  were  included.  Two  replications  of  each  test  were 
completed,  except  that  three  replications  were  run  for  the  three  seconds,  75  kW/m2  exposures. 

All  garments  were  one-piece  coveralls  of  a similar  design.  Garments  were  produced  by  three 
different  manufacturers  and  were  therefore  cut  somewhat  differently.  Garment  A was  size  42R 
and  fit  somewhat  snugly  on  the  mannequin.  Garment  B was  size  44  R and  fit  the  mannequin 
loosely.  Garments  C,  D,  E and  F were  size  42R  and  fit  the  mannequin  well.  All  garments  were 
washed  once  in  liquid  Tide  to  remove  any  residual  mill  finishes  before  testing.  Laundry  conditions 
included  50°C  wash  for  10  minutes,  cold  rinse  and  either  drying  in  the  dryer  for  20-30  minutes 
on  the  permanent  press  cycle  (for  garments  D,  E and  F)  or  hanging  to  dry  (garments  A,  B and  C). 

Characteristics  of  the  garment  fabrics,  determined  using  extra  fabric  samples  supplied  by  the 
garment  manufacturer,  are  outlined  in  Table  3 with  the  standard  test  procedures  noted.  Fabric 
mass  was  determined  prior  to  laundering;  fabric  count  and  thickness  were  determined  after 
laundering.  Thermal  protective  performance  (TPP)  was  measured  after  laundering  following 
ASTM  D4108,  as  modified  in  Par.  6.1  of  CAN/CGSB-155.1-M88.  This  procedure  uses  an  open 
flame  and  no  spacer  between  the  fabric  and  the  calorimeter.  The  fabric  is  held  in  place  in  the 
specimen  holder  by  a set  of  pins. 


TABLE  3 Fabric  Characterization  and  Thermal  Protective  Performance  Ratings  of  Fabrics  Used  in  Garments  Tested 
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RESULTS 


The  percent  of  mannequin  surface  reaching  the  2nd  or  3rd-degree  bum  criteria  is  reported 
for  each  replication  of  each  test  in  Table  4,  demonstrating  considerable  consistency  between 
replications  for  any  exposure  condition.  The  mean  value  for  each  garment  type  for  each  heat 
flux/exposure  time  combination  is  given  in  Fig.  7.  All  values  include  7%  for  the  unprotected 
Lead  of  the  mannequin.  Since  there  are  no  sensors  in  the  hands  or  feet,  however,  the  maximum 
possible  area  reaching  the  bum  criteria  is  86%. 

Typical  computer  print-outs  of  results  for  garments  B and  E,  exposed  at  75  kW/m2  (1.8 
cal/cm2*s)  for  both  three  and  four  seconds  are  shown  in  Figs.  8 to  11.  The  figures  differentiate 
between  the  2nd  and  3rd-degree  criteria.  A typical  plot  of  % burn  with  time  is  shown  in  Fig.  12. 

The  effect  of  varying  heat  flux  and  exposure  time  differs  with  garment  type.  While  the 
heavier  fabrics  with  the  highest  TPP  values  seem  to  protect  very  well  at  the  lowest  and  shortest 
exposures,  the  lightest  garment  with  the  lowest  TPP  value  (Garment  B)  gave  the  best  protection 
from  the  most  severe  exposure  (4  seconds,  75  kW/m2).  Visual  inspection  of  the  garments  after 
exposure  indicate  that  while  Garments  A,  B and  C kept  their  fabric  integrity  throughout  the  most 
severe  exposure.  Garments  D,  E and  F did  not.  Garments  D and  E experienced  considerable  after 
flame  during  this  exposure,  and  when  the  researchers  attempted  to  remove  these  garments  from 
the  mannequin,  they  fell  apart.  The  fabric  integrity  of  Garment  F was  affected  less  than  that  of 
garments  D and  E.  Following  exposure,  a crystalline  material  was  evident  on  the  inside  of 
Garment  F and  on  the  mannequin  surface. 

These  results  suggest  that  the  choice  of  exposure  conditions  greatly  affects  the  relative 
results  of  garment  testing  on  an  instrumented  mannequin.  Thus,  whenever  possible,  exposure 
conditions  should  be  selected  to  most  closely  simulate  the  hazardous  environment  from  which 
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protection  is  sought  The  results  also  draw  attention  to  a limitation  of  the  TPP  tests  in  common 
use.  Further  investigation  of  fabric  integrity  when  exposed  to  high  heat  flux  for  exposure  over 
three  seconds  is  recommended. 


TABLE  4 Percent  of  Mannequin  Surface  Reaching  2nd  or  3rd-Degree  Burn  Criteria21 


Garment 

Heat  Flux 

67  kW/m2 

75  kW/m2 

84  kW/m2 

TPP 

(Rep) 

Exp.  Time 

3 sec 

4 sec 

3 sec 

4 sec 

3 sec 

Mean  (St’d  dev.) 

A (1) 

15.3 

46.8 

29.0 

63.3 

45.2 

5.2 

(2) 

16.8 

56.6 

37.4 

66.1 

35.8 

(±0.3) 

(3) 

33.6 

B (1) 

10.4 

30.5 

12.3 

47.3 

17.5 

4.3 

(2) 

10.4 

24.7 

13.1 

47.5 

13.9 

(±0.4) 

(3) 

11.7 

C(l) 

8.4 

28.7 

17.5 

55.1 

28.7 

4.9 

(2) 

8.3 

30.4 

11.7 

51.9 

26.6 

(±0.3) 

(3) 

18.2 

D(l) 

8.5 

34.5 

9.8 

70.8 

10.2 

7.3 

(2) 

7.9 

26.4 

11.3 

75.5 

10.3 

(±0.5) 

(3) 

9.9 

E(l) 

8.5 

37.6 

9.8 

80.4 

13.3 

7.3 

(2) 

8.8 

24.4 

10.7 

81.5 

9.7 

(±0.4) 

(3) 

8.8 

9.1 

F (1) 

14.4 

32.1 

14.7 

58.3 

15.7 

6.4 

(2) 

11.6 

32.7 

17.3 

65.7 

13.9 

(±0.2) 

(3) 

13.4 

aMaximum  possible  % is  86,  since  hands  and  feet  do  not  have  sensors;  figures  include  1% 
unprotected  mannequin  in  the  head  area: 


Garments 


A:  210  g (6.2  oz)  plain  weave  Nomex®  HI  (aramid),  size  42R  (close-fitting) 
B:  187  g (5.5  oz)  plain  weave  60/40  Kevlar®  (aramid)/Pbi®,  size  44 
C:  249  g (7.4  oz)  rib  weave  Nomex®  IH  (aramid),  size  42R 
D:  307  g (9.1  oz)  satin  weave  Proban®  FR-7A  cotton,  size  42R 
E:  334  g (9.9  oz)  satin  weave  Indura®  Proban®  cotton,  size  42R 
F:  305  g (9.0  oz)  twill  weave  Vulkan®  (FR  viscose/aramid),  size  42R 
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Fig,  7 Mean  percent  mannequin  surface  reaching  2nd  or  3rd  degree  bum  criteria 

(Note:  Maximum  possible  % is  86,  as  hands  and  feet  do  not  have  sensors;  all  figures 
include  7%  unprotected  mannequin  in  the  head  area.) 
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Garment  B,  75  kW/m2,  4 second  exposure 
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Fig.  1 1 Calculated  location  of  skin  burning  on  mannequin 
Garment  E,  75  kW/m2,  4 second  exposure 
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Fig.  12  Calculated  cumulative  skin  burning  for  a typical  4 second  exposure 
with  protective  clothing 
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D.  SUMMARY  AND  CONCLUSIONS 

A full-sized  instrumented  mannequin  has  been  developed  which  can  be  used  to  evaluate  the 
thermal  protective  qualities  of  garments  in  flash  fires  or  fire  balls  of  varying  intensity  and 
duration.  Garment  testing  has  shown  that,  in  general,  increasing  the  fabric  weight  and/or 
thickness  for  a particular  fiber  or  fiber/finish  combination  reduces  the  extent  of  skin  burning. 
However,  fabrics  of  some  polymeric  materials  (e.g.  Pbi®)  offer  more  protection  than  others,  so 
that  a general  correlation  between  thickness/weight  and  protection  does  not  hold,  especially  at 
higher  heat  fluxes  and  durations.  TPP  values  correlate  better  with  the  predicted  amount  of  skin 
burning  on  the  mannequin  at  low  heat  flux/duration  than  at  higher  heat  flux/duration  exposures. 

The  results  of  this  research  indicate  that  for  garment  testing,  exposure  conditions  should 
be  selected  to  simulate  actual  flash  fire/fireball  hazards.  The  effect  of  fabric  integrity  on 
protection  at  high  exposures  (both  heat  flux  and  duration)  needs  further  investigation. 
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APPENDIX 


A COMPUTER  CODE  FOR  ESTIMATING  THERMAL  INJURY  OF  HUMAN  SKIN 

This  program  was  written  in  QuickBasic  4.5™  for  using  data  from  an  instrumented 
mannequin  to  estimate  the  nature  and  amount  of  skin  burning  when  the  mannequin  is  exposed 
to  flash  fires.  In  general,  the  program  uses  Duhamel’s  Integral  Method  to  calculate  the  surface 
heat  flux  from  the  measured  surface  temperatures  of  the  sensors.  On  the  basis  of  the  calculated 
surface  heat  flux,  a multi-layer  skin  model  is  solved  for  the  skin  temperatures  using  the  Crank- 
Nicolson  Implicit  method.  From  the  skin  temperature  distribution,  the  degree  of  thermal  injury 
and  the  depth  of  damage  for  each  location  of  the  body  are  determined  by  a first  order  reaction 
rate  process  equation  coupled  with  a bum  criterion  from  the  literature.  The  percentage  of  the 
total  area  burned  is  then  estimated  based  on  the  number  of  sensors  which  indicate  2nd  degree 
bum  or  above.  The  program  also  draws  the  burned  area  of  the  mannequin  at  each  second  of 
realtime  up  to  any  chosen  time  period  (max.  2 minutes)  and  plots  the  percent  cumulative  bum 
versus  time.  In  this  program  integrations  of  all  time  functions  are  performed  using  the 
trapezoidal  rule  and  the  thermal  response  of  the  skin  is  assumed  to  be  independent  of  site.  The 
above  procedures  are  described  in  Fig.  Al. 


Fig.  A1  The  Flow  Diagram  of  the  Program 
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To  lead 
Data  File  5 
again  for  the 
next  sensor 

A 


From  Sldn  Model 


Born  Model 

From  the  «ldn  temperature  distribution*  this  model 

1.  Calculates  the  damage  function  Q at  different 
depths  for  each  location  by  using  one  of  the 
available  bum  criteria  horn  the  literature,, 

2.  Converts  Q into  1st,  2nd  or  3rd  degree  bum 

and  determines  the  depth  of  damage  and  the 
times  to  2nd  and  3rd  degree  bums. 


Determination  of  Percent  Burned  Area 

The  bum  prediction  at  the  end  of  the 
period  of  study  for  each  sensor  is  added 
together  to  find  the  percentage  of  area 
which  has  1st,  2nd  or  3rd  degree  bum. 


v 


Read  Data  FTT  B $ 

for  the  coordinates  of  points  contouring  the 
mannequin  and  the  areas  of  the  sensors. 




Draw  Man 

Based  on  the  times  to  2nd  and  3rd  degree  bums 
obtained  from  the  Bum  Model  for  each  of  the 
sensor  locations,  the  % cumulative  2nd  and  3rd 
degree  bums  with  respect  to  time  are  calculated, 
and  the  areas  of  the  mannequin  which  have  at 
least  2nd  degree  bum  are  drawn  on  the  screen 
for  each  sec.  interval. 


_Y_ 

Plot  die  % cumulative  bums  versus 
time  on  the  screen. 


Fig.  A1  (continued) 
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V 


Output  FT1VR  7 

Contains  the  time  - heat  flux  - sensor  surface  temperature  - 
skin  temperature  distribution  data,  and  die  damage  function 
Q at  the  basal  layer  and  at  the  depth  of  damage  for  any 
chosen  sensor  locations  (optional). 
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V 


Output  E2LE  $ 

Contains  the  summary  of  sensor  response  and  the  mannequin 
bum  prediction.  The  former  includes  the  times  to  2nd  and 
3rd  degree  bums,  the  total  accumulated  heats  at  2nd  degree 
bum  and  at  the  end  of  the  period  of  study,  the  depth  of 
damage  and  the  predicted  degree  bum  for  each  sensor  The 
latter  includes  the  predicted  percent  areas  of  1st,  2nd  and  3rd 
degree  bums  for  the  whole  mannequin,  and  the  sums  of  the 
total  accumulated  heats  for  all  of  the  sensors. 


Fig.  A1  (continued) 
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The  input  data  required  for  this  program  are  stored  in  5 data  files.  File  1 contains  mainly 
the  thermal  inertia  (kpC)  of  the  sensors.  The  data  in  File  1 are  formatted  as  follows: 


NSEN 

NOSj 

KPCj 

SLOPEj 

INTERCEPT! 

nos2 

kpc2 

slope2 

INTERCEPT 

N0S(NSEN) 

^^(NSEN) 

slope(NSEN) 

INTERCEPT  (nSen) 

Where: 

NSEN  = Total  number  of  sensors 

NOSj  = A sensor  number  which  identifies  the  location  of  the  sensor  on  the  mannequin. 

The  subscript  i indicates  the  i-th  sensor  on  the  file. 

KPCa  = The  thermal  inertia  of  the  sensor  identified  by  NOS;,  cal2/cm4*°C2*sec. 

SLOPE!  = The  slope  of  the  calibration  line  for  the  sensor  identified  by  NOSi? 

cm2*°C*seca5/caL 

INTERCEPT!  = The  intercept  of  the  calibration  line  for  the  sensor  identified  by  NOS^ 
°C-cm2-sec/cal. 

File  2 contains  the  percent  area  each  sensor  occupies  (PCT)  which  are  formatted  as: 

NOSj  PCTj  NOS2  PCT2 PCT(nSEn) 

File  3 is  required  only  for  drawing  the  mannequin  and  its  burned  area.  Thus  File  3 consists  of 
the  coordinates  which  contour  the  mannequin  and  the  areas  of  the  sensors.  These  data  are 


formatted  as: 
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Where: 

N1 

FXi 

FT, 

N2 

BXi 

BY, 

N3 

FXMANi 

FYMANX 

N4 

BXMANi 

BYMANj 


N1 

FXj 

fx2 

FX^d 

FY, 

fy2 

FY^j) 

N2 

BXj 

BX2 

bx(N2) 

BYj 

by2 

BY^) 

N3 

FXMANj 

fxman2 

....  fxman(N3) 

FYMANj 

fyman2 

FXMAN^) 

N4 

BXMANj 

bxman2 

....  BXMAN^ 

BYMANj 

byman2 

....  BYMAN^ 

= No.  of  points  required  to  outline  the  front  view  of  the  mannequin 
= X-coordinate  for  the  front  view  of  the  outline  of  the  mannequin 

= Y-coordinate  for  the  front  view  of  the  outline  of  the  mannequin 

= No.  of  points  required  to  outline  the  back  view  of  the  mannequin 
= X-coordinate  for  the  back  view  of  the  outline  of  the  mannequin 

= Y-coordinate  for  the  back  view  of  the  outline  of  the  mannequin 

= No.  of  points  required  to  outline  the  front  view  of  the  areas  of  the  sensors 
= X-coordinate  for  the  front  view  of  the  areas  of  the  sensors 

= Y-coordinate  for  the  front  view  of  the  areas  of  the  sensors 

= No.  of  points  required  to  outline  the  back  view  of  the  areas  of  the  sensors 
= X-coordinate  for  the  back  view  of  the  areas  of  the  sensors 

= Y-coordinate  for  the  back  view  of  the  areas  of  the  sensor 


and,  subscript  i = the  i-th  coordinate  point. 
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In  File  4,  all  the  thermal  properties  of  human  skin  and  the  control  parameters  are 


included. 


Where: 
Mode  1 

Mode  2 


DT 

NST 

NG1 

NG2 

NG3 

NS1 

NS2 

NS3 

R1 

R2 

R3 

TI 


The  format  of  the  data  in  this  file  is  shown  as  follows: 


Mode  1 

Mode  2 

DT 

NG1 

NG2  NG3 

NS1 

R1 

R2  R3 

TI 

TB  CPB 

H 

THK(l) 

THK(2) 

THK(3) 

CP(1) 

CP(2) 

CP(3) 

COND(l) 

C0ND(2) 

C0ND(3) 

CKD 

G(2) 

G(3) 

GAMA(l) 

GAMA(2) 

NST 

NS2 


NS3 


= A mode  number  for  selecting  the  type  of  skin 

1.  for  opaque  skin  (ignore  penetrating  radiant  heat) 

2.  for  diathermanous  skin  (include  penetrating  radiant  heat) 

= A mode  number  for  selecting  the  bum  criteria 

1.  for  Henriques’  criteria 

2.  for  Stoll’s  criteria 

3.  For  M.I.T.’s  criteria 

4.  combined  criteria  of  Stoll’s  for  epidermal  burn  and  Takata’s  for  dermal  layer 
= time  step,  sec. 

= printing  steps  for  time 

= No.  of  grid  in  epidermis 

= No.  of  grid  in  dermis 

= No.  of  grid  in  subcutaneous  tissue 

= printing  steps  for  NG1 

= printing  steps  for  NG2 

= printing  steps  for  NG3 

= axial  stretching  ratio  for  epidermis 

= axial  stretching  ratio  for  dermis 

= axial  stretching  ratio  for  subcutaneous  tissue 

= initial  temperature  at  skin  surface,  °C 
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TB  = blood  temperature,  °C 

CPB  = volumetric  heat  capacity  of  blood,  cal/cm3  °°C 
H = overall  surface  heat  loss  coefficient,  cal/cm2  *sec-°C 

THK(k)  = thickness  of  the  k-th  layer  of  skin,  cm 
Cp(k)  = volumetric  heat  capacity  of  the  k-th  layer  of  skin,  cal/cm3  *°C 
COND(k)  = conductivity  of  the  k-th  layer  of  skin,  cal/cm -°C -sec 
G(k)  = blood  perfusion  rate  for  skin  layer  k,  cm3/sec  per  cm3  of  tissue 
GAMA(l)  =absorption  coefficient  of  epidermis,  cm"1 
GAMA(2)  =absorption  coefficient  of  dermis  and  subcutaneous  tissue,  cm"1 
and,  k = 1 for  epidermis 

= 2 for  dermis 

= 3 for  subcutaneous  tissue 

The  last  data  file,  File  5,  stores  all  the  measured  surface  temperatures  of  the  sensors  at 
each  time  interval  in  the  following  format: 

NOS(l) 

TSS(l)  TSS(2)  ...  999.9 

NOS  (2) 

TSS(l)  TSS(2)  ...  999.9 


NOS(NSEN) 

TSS(l)  TSS(2)  ...  999.9 

Where: 

NOS(M)  = The  sensor  number  of  the  M-th  sensor  the  file  stores  its  temperatures 
TSS(J)  = The  measured  sensor  surface  temperature  at  the  J-th  point  in  time,  °C 

The  number  999.9  at  the  end  of  each  temperature  line  is  used  to  indicate  the  end  of  the 


data  for  one  sensor. 
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In  addition  to  the  input  data  stored  in  the  files  as  mentioned  previously,  the  program  also 
requires  some  other  information.  This  information  is  input  directly  from  the  keyboard  and 
include  the  following: 


COVER 

= Type  of  clothing  system 

WEIGHT 

= Fabric  mass,  g/m2 

FABTHK 

= Fabric  thickness,  mm 

XPS 

= Exposure  heat  flux,  cal/cm2<sec 

XPSTME 

= Exposure  time,  sec 

NSENP 

= Number  of  sensors  selected  for  printing 

NOSP 

= Selected  sensor  numbers  for  printing 

TA 

= Ambient  temperature,  °C 

TMC 

= Calculation  time,  sec 

In  the  following  sections,  the  sensor  model,  the  skin  model  and  the  burn  model  will  be 


discussed  in  detail. 
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SENSOR  MODEL 

This  model  assumes  that  the  sensor  (the  skin  simulant)  is  a semi- infinite  solid  with 
constant  thermal  properties  and  the  heat  transfer  phenomena  obey  the  one-dimensional  Fourier 
equation: 


C = specific  heat  capacity,  cal/g*°C 

These  are  valid  assumptions  for  the  short  time  periods  involved  with  flash  fires. 

Given  the  measured  surface  temperature,  this  model  calculates  the  surface  heat  flux  of 
the  sensor  using  the  Duhamel’s  Integral  Method. 

The  Duhamel’s  theorem  states  that  if  u^x,!)  is  the  response  of  a linear  homogeneous 
system  (initially  at  zero)  to  a single  unit  step  input,  the  response  uF(x,t)  of  the  system  to  the  input 
F(t)  is  given  by  the  following  equation: 


a 2t 
d? 


(Al) 


where  T =temperature,  °C 


t =time,  sec 


x =distance  below  the  surface,  cm 


a = k/pC  - thermal  diffusivity,  cm2/sec 


k ^thermal  conductivity,  cal/cm*sec*°C 
p = density,  g/cm3 


(A2) 


o 


In  order  to  calculate  the  surface  heat  flux  based  on  the  measured  surface  temperature  using  the 
Duhamel’s  theorem,  one  has  to  know  how  the  surface  heat  flux  responds  to  a single  unit  step 
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increase  in  surface  temperature.  The  classical  solution  of  Eq.  (Al)  for  the  temperature  response 
of  a semi-infinite  solid  initially  at  Tj  and  subjected  to  a step  change  in  surface  temperature  Ts  at 
zero  time  has  been  given  in  the  literature  [A6]: 


T-T: 


‘ = 1-JL  (W“  e-y‘dy 

& Jo 


T -T 

S l 


(A3) 


Thus  the  surface  heat  flux  q is: 


9(0  = -k  (Ts-T;)  JL  Jf2^  e'y2dy 


= 31  (Ts-Tt) 

\TZ 


2k 


_d(x/2<JaT ) 
2/4w  1 


jxf2]/ai  e-y2dy 


dixllJaT ) 


dx 


2\[aT 


*=0 


7 zat 


(A4) 


For  a single  unit  step  increase  in  surface  temperature,  Ts-T1=  1.  Thus,  the  heat  flux  response 
(initially  at  zero)  becomes: 


<7(0  = 


\j%at 


(A5) 


Substituting  Eq.  (A5)  into  Eq.  (A2),  the  corresponding  heat  flux  for  a variable  surface 
temperature  is: 
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Q(f ) = J W-r,.) 


fix 


i (*-x)3/2  i (?-x)3/2 


2(0  = 


kpC 

n 


Ts(t)~Ts(x) 

(?-x)3/2 


fix  + 


Ts(t)~Tt 

t(U2) 


(A6) 


The  integral  in  Eq.  (A6)  can  be  evaluated  easily  by  standard  numerical  integration 
methods  such  as  the  trapezoidal  rule. 

One  criterion  for  choosing  the  skin  simulant  is  that  its  thermal  inertia,  kpC,  should  match 
that  of  the  skin  [A3].  To  find  the  thermal  inertia,  the  surface  of  the  simulant  was  irradiated  with 
a known  constant  heat  flux,  and  the  surface  temperature  rise  (AT)  with  time  was  measured.  The 
thermal  inertia  can  then  be  obtained  from  Eq.  (A7),  given  in  Reference  [A6],  describing  the 
temperature  response  of  a semi-infinite  solid  initially  at  Tx  to  a constant  heat  flux  q at  t = 0: 


T(x,t)-Ti  = « 
k 


_ X2 

f \ 

2 

e 4ctr  - x etfc 

X 

\ 

7t 

2\[aT  J_ 

At  x = 0,  Eq.  (A7)  becomes: 


(A7) 


AT  2 
Q \lkpCiz 

Thus,  by  plotting  AT/q  vs.  Vt,  kpC  can  be  calculated  from  the  slope  2/VkpC7t. 

The  skin  simulant  is  an  inorganic  mixture  of  calcium,  aluminum,  silicate  with  asbestos 
fibres  and  a binder,  which  is  called  "Colorceran".  Each  sensor  which  is  made  up  of  this  material, 
is  about  19  mm  in  diameter  and  32  mm  long.  To  measure  the  surface  temperature  a 30  gauge 


s/T 


(A8) 
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Cu-Co  thermal  couple  wire  rolled  flat  to  about  0. 1 mm  thick  was  held  on  the  surface  with  an 
epoxy-phenolic  adhesive.  Figure  A2  shows  the  plot  of  Eq.  (A8)  for  a test  sensor.  The  thermal 
inertia  calculated  from  the  slope  is  about  12.53  x 10"4  cal2/cm4*°C2*sec.  The  density  of  this 
material  was  measured  to  be  1.877  g/cm3.  The  heat  conductivity  obtained  from  an  independent 
measurement  [A7]  was  approximately  0.97  W/mK  (23.18  x 10"4  cal/cm-sec *°C).  Once  these 
thermal  properties  are  known,  the  specific  heat  capacity  and  diffusivity  can  easily  be  calculated. 
These  are  summarized  in  Table  A1  in  which  the  thermal  properties  of  human  skin  [3,8]  are  also 
included. 

TABLE  A1  Thermal  Properties  of  Human  Skin  and  Skin  Simulant 


Properties 

Human  Skin 

Skin 

Simulant 

Epidermis 

Dermis 

k(cal/cm*sec*°C) 

6.i  x nr4 

12.5  x 10-4 

23.18  x nr4 

p(g/cm3) 

1.2 

1.2 

1.877 

C(cal/g-°C) 

0.86 

0.77 

0.288 

a(cm2/sec) 

5.91  x 10-4 

13.53  x 10-4 

42.88  x 10-4 

kpC(cal2/cm4-°C2-sec) 

6.30  x 10-4 

11.55  x 10-4 

12.53  x 10-4 

As  seen  in  Table  Al,  the  thermal  inertia  of  the  skin  simulant  is  within  10%  of  that  of  the 
dermis  which  occupies  about  96%  of  the  thickness  of  the  skin.  Therefore,  "Colorceran"  is 
appropriate  for  using  as  a skin  simulant. 

There  were  1 10  of  these  skin  simulant  sensors  installed  on  the  mannequin.  Each  of  these 
sensors  was  calibrated  using  the  same  method.  The  variation  of  the  thermal  inertia  for  these 
sensors  are  shown  in  Fig.  A3.  As  indicated,  the  mean  thermal  inertia  is  11.71  x 10"4 
cal2/cm4*°C2*sec  with  a standard  deviation  of  1.11  x 10-4. 


AT/Q  (°C  cm2  sec/cal) 
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Fig.  A2  Calibration  result  of  the  thermal  inertia  of  a skin  simulant  sensor 


Number  of  Sensors 
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50 


Mean  kpc  = 1 1.71  x 10"4  cal2/cm4  °C2sec 
Standard  deviation  =1.11  x TO"4  cal2/cm4  °C?sec 
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Fig.  A3  Variation  of  the  thermal  inertia  of  the  110  skin  simulants 
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Due  to  the  surface  condition  of  the  sensors  and  the  adhesion  situation  between  the 
thermal  couples  and  the  surfaces,  the  plots  of  Eq.  (A8)  usually  do  not  intercept  exactly  at  the 
origin.  The  mean  Y-intercept  of  the  110  sensors  was  found  to  be  -1.5326  °Ct;m2‘sec/cal  with 
a standard  deviation  of  1.989.  The  negative  value  of  the  mean  Y-intercept  indicates  a time  lag 
in  initial  response  for  the  majority  of  the  sensors  and  the  time  lag  is  equal  to  (Y-intercept/slope)2. 

Since  the  Duhamel’s  Integral  Method  is  only  applicable  to  a linear  homogeneous  system 
of  which  the  response  to  an  input  function  is  initially  at  zero  (i.e.  for  Y-intercept  equal  zero),  the 
calculation  of  the  heat  flux  using  this  method  by  neglecting  the  Y-intercept  would  lead,  to  a 
computational  error.  Thus  a simple  method  was  derived  to  correct  the  prediction  of  heat  flux 
using  the  Duhamel’s  Integral  Method  without  consideration  of  the  Y-intercept. 

Recall  Eq.  (A8).  That  is: 

= slope  yJT  (A8) 

Q 

where  slope  - 2/VkpC7t.  With  intercept  included,  Eq.  (A8)  becomes: 


AT  , r~ 

— - slope  \Jt  + intercept 

Q 

To  predict  the  heat  flux  (qj)  without  consideration  of  the  intercept: 


(0  = 


AT 


slope  \[t~ 

To  predict  the  heat  flux  (q2)  with  consideration  of  the  intercept: 


(A8a) 


<?2(0 


AT 


slope  yfT  + intercept 


(A8b) 
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To  predict  q2  from  qv  combine  Eqs,  (A8a)  and  (A8b): 


_ slope  \jT 

slope  \[T  + intercept 

so  that  q2(t)  = CF  x qx(t)  where 


1 

^ + intercept  1 
slope  Jf 


1 + intercept  1 (A9) 

slope 

As  indicated  in  Eq.  (A9),  when  time  increases,  CF  increases  if  the  intercept  is  positive,  or 
decreases  if  the  intercept  is  negative.  When  time  tends  to  infinity,  CF  tends  to  one.  Equation 
(A9)  also  indicates  that: 

When  t > 100  (intercept/slope)2, 


* >0.91  and  < 1 for  positive  intercept 
> 1 and  <1.11  for  negative  intercept 

For  the  test  results  shown  in  Fig.  A3,  the  mean  slope  and  intercept  are  about  33.0 
cm2*°C*  sec,/2/cal  and  -1.536°C*  cm2*  sec/cal  respectively.  If  a sensor  has  the  mean  slope  and 


intercept,  the  computational  error  of  the  heat  flux  predicted  by  neglecting  the  intercept  would 
drop  to  about  10%  in  0.22  sec.  Therefore,  Eq.  (A9)  is  significant  only  at  the  beginning  of  time. 
However,  when  t < (intercept/slope)2,  Eq.  (A9)  becomes: 


CF 


< >0 
.<  o 


and 


< 0.5  for  positive  intercept 
for  negative  intercept 


Practically  CF  can  never  be  less  than  zero.  Therefore,  for  the  case  of  negative  intercept,  Eq. 
(A9)  can  only  be  used  for  t > (intercept/slope)2,  i.e.  for  time  greater  than  the  time  lag.  The 
prediction  of  heat  flux  for  each  sensor  was  corrected  using  its  unique  correction  factor,  CF. 
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SKIN  MODEL 

The  model  used  was  developed  by  Mehta  and  Wong  from  M.I.T.  [A3].  The  skin  was 
modelled  as  consisting  of  three  layers:  the  epidermis,  the  dermis  and  the  subcutaneous  tissue, 
each  of  uniform  thickness  as  shown  in  Fig.  A4a.  The  arteries  and  veins  are  assumed  to  be 
straight  to  the  dermis-epidermis  interface.  The  deepest  subcutaneous  layer  is  considered  to  be 
at  the  constant  normal  body  temperature  (37 °C).  The  blood  in  the  artery  is  assumed  to  enter  at 
this  temperature.  As  the  blood  travels  outward  from  the  isothermal  core,  part  of  it  will  leave  the 
arteries  through  the  capillaries  to  the  veins  and  finally  back  to  the  isothermal  core.  In  this  model 
each  of  the  three  layers  is  assumed  to  have  constant  thermal  properties.  The  blood  perfusion  rate 
is  assumed  to  be  constant  within  the  lower  two  layers,  but  is  zero  within  the  epidermis.  The  skin 
can  be  considered  opaque  or  diathermanous  depending  on  the  source  of  thermal  radiation. 
Opaque  skin  ignores  the  effect  of  penetrating  radiant  heat  flow,  while  the  diathermanous  skin 
model  accounts  for  the  penetration  effects  of  thermal  radiation. 

In  developing  the  energy  balance  equations  for  the  three  layers,  Mehta  and  Wong  further 
assumed  that: 

1.  The  heat  transfer  is  a one-dimensional  transient  flow  through  the  skin  layers. 

2.  The  outlet  temperature  of  the  blood  from  the  capillaries  to  the  veins  equilibrates  with  the 
local  tissue  temperature  because  of  the  large  surface-to-volume  ratio  and  relatively  long 
time  of  residence  of  blood  in  the  capillaries. 

3.  Heat  transfer  between  the  larger  blood  vessels  (arteries  and  veins)  and  the  tissue  is 
negligible.  The  temperature  of  the  arterial  blood,  unlike  the  blood  temperature  in  the 
veins,  is  independent  of  x and  equal  to  the  normal  body  temperature. 

4.  Transmission  of  radiation  through  diathermanous  skin  obeys  the  Beer-Lamber  law  with 
one  absorption  coefficient  applying  for  the  epidermal  layer,  and  another  for  the  deeper 
tissue.  Absorption  of  radiation  by  the  blood  is  negligible. 


Heat  Source 
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x 


Fig.  A4a  The  skin  model 


50 


Arteries 


Fig.  A4b 


Capillaries 


A volume  element  of  skin  indicating  the  blood  flow  (solid  arrows)  and  heat 
transfer  (dotted  arrows) 
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On  the  basis  of  the  above  assumptions,  the  transient  flow  of  heat  through  a volume  tissue 
element  as  shown  in  Fig.  A4b  can  be  described  by  the  following  equation: 

Cpil=  k^L  - (Cp)bG(T-Tb)  + qf  (A10) 

ot  dx2 

The  left  hand  side  of  Eq.  (A  10)  represents  the  rate  of  increase  of  energy  content  of  the 
tissue  in  which  Cp  is  the  volumetric  heat  capacity  in  cal/cm3*°C  and  T is  the  tissue  temperature. 
The  first  term  on  the  right  side  represents  the  conductive  heat  gain  within  the  tissue,  where  k is 
the  thermal  conductivity  in  cal/cm*  sec*°C.  The  second  term  on  the  right  is  the  heat  transfer 
from  the  tissue  to  the  capillaries  in  which  (Cp)b  denotes  the  volumetric  heat  capacity  of  the 
blood,  G is  the  blood  perfusion  rate  in  cc/sec  per  cc  of  tissue  volume  and  Tb  is  the  perfusing 
blood  temperature  (i.e.  the  arterial  blood  temperature,  but  assumed  to  be  at  37 °C,  the  internal 
body  temperature).  The  last  term  represents  the  absorption  of  penetrating  radiation  by  the  tissue 
in  which  q is  the  absorbed  radiation  at  the  surface  in  cal/sec*  cm2  and  f denotes  the  energy 
attenuation  functions  which  are  defined  as  follows: 


/ = /M  = Yi  e 0 < * < 5j 


(All) 


= *>6, 

where  y1  = absorption  coefficient  for  the  epidermal  layer,  cm'1 

y2  = absorption  coefficient  for  the  dermal  and  subcutaneous  layers,  cm'1 

= thickness  of  epidermis,  cm. 

Equation  (A  10)  can  be  solved  for  the  three  layers.  For  opaque  skin,  the  last  term  in 
Eq.  (A  10)  is  omitted.  In  the  epidermal  layer,  the  blood  perfusion  G is  zero  because  of  its 
nonvascularity.  The  initial  and  boundary  conditions  to  be  satisfied  are: 
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(i)  t = 0,  T = T0(x)  = an  initial  temperature  profile 

(ii)  t > 0,  x = 0, 

-kx  = q(t)  - h(T-T,J  for  opaque  skin 
dx 


rlT 

or  kj  __  = hOT-T^  for  diathermanous  skin 
dx 

where  kj  = thermal  conductivity  of  the  epidermis,  cal/sec-  cm-°C 
q(t)=  radiation  heat  flux  at  the  surface,  cal/cm2*  sec 
h = overall  surface  heat  loss  coefficient,  cal/cm2'  sec*°C 
= ambient  temperature,  °C 

(iii)  t > 0,  x = 83,  T = Tb 

where  53  = distance  of  the  deepest  subcutaneous  layer  from  the  surface,  cm 
Tb  = normal  body  temperature,  °C 

(iv)  Continuity  of  temperature  and  heat  flux  at  the  two  interfaces  at  all  times. 

To  solve  Eq.  (A  10)  numerically,  a non-uniform  grid  system  was  devised  with  fine  grid 
near  the  skin  surface  and  the  two  interfaces,  which  was  stretched  to  a coarse  size  in  the  regions 
where  the  temperature  gradients  are  smaller.  This  grid  system  is  shown  in  Fig.  A5  and  described 
in  the  following  equations. 


for  i = 1,2,3,...,N7-1 

Ajj 

for  i = N1 

A n (i-Nl-1) 

^N1 +1*'2 

for  i = N1  +1,  N1  +2 +A2  -1 

A* A/i+1 

for  i = NJ  +N2 

A p(i-Nl-N2-l) 

^Nl +N2 +lK3 

for  i = N1  +N2  +N2  +N3 

(A  12) 


53 


Fig.  A5  A nonlinear  finite  difference  grid  for  the  skin  model 
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where 


Ax ! = 


Ax 


N1+ 1 


N7-1 

1 + Y, 

i’=l 

^2”^1 
M+A/2-1 
1+  £ 

/=/V7  +1 


Ax 


83  82 


N7+/V2+1 


NJ  +N2  +N3 


E ^3 

i=Nl*N2+ 1 


(HW-N2-1) 


N1  = number  of  grids  in  epidermis 

N2  = number  of  grids  in  dermis 

N3  = number  of  grids  in  subcutaneous  tissue 

Rj  = axial  stretching  ratio  for  epidermis 

R2  = axial  stretching  ratio  for  dermis 

R3  = axial  stretching  ratio  for  subcutaneous  tissue 

Also,  to  solve  Eq.  (A  10)  numerically,  it  is  more  convenient  to  express  Eq.(AlO)  in  a 
dimensionless  form  by  defining: 


e = (T  - Tb)/Tb 

(A13) 

> 

II 

X 

(A  14) 

0 

> 

II 

(A  15) 

<J>  = (to/CpTb)qf 

(A  16) 

where  t0  = exposure  time. 
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Thus,  Eq.  (A  10)  becomes: 


30 

a7 


o a2e  , 

(3__  - X0  + 

ax2 


where  p 


_L  !l 
Cp  8? 


X. 


(A  17) 


Based  on  the  concept  of  the  Crank-Nicholson  Implicit  Method,  each  term  in  Eq.  (A  17) 
is  written  in  a finite  difference  form  at  the  point  (i,  j+Vi)  in  the  time-distance  coordinate  plane 
as  follows: 


i-l,j  i,j  i+l»j 


Axi-1 

AX. 

l 

— 2 
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•1 
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90  | = e;,;*i  -ey 

W‘-^A  " Ax 


92e 


ft 


dx2  'y+*  " ax,+axm 

Mlu*  - ~ [VeyJ 


®i,j 


AX, 


e*V"0i-l  J 
AX,_x 


®/+i,y+i  ®i,y+i 


AX; 


®/,y+i  ®i-i,y+i 
AX^ 


O-  • = J_  [d>.  +<i> . J 
ij+^  ^7  L JJ+1J 


where  the  subscripts  i and  j denote  the  points  in  coordinates  X (or  x)  and  x (or  t)  respectively, 
and  start  at  one. 

The  dimensionless  finite  different  form  of  the  initial  and  boundary  conditions  can  be 

expressed  as: 

(i)  0U  = (Tu  - Tb)/Tb  i = 1,  2,  ....  N1+N2+N3+1  (A18) 

where  T,  x = T0(xj)  = the  initial  temperature  profile 


(ii)  6J±Sl  = H[0,  j - 6J  - <j>;,  j = 2,  3,  4, ... 
2AX,  1 1 


where  H = hbj/kj 


5i  q. 

(t>j  = J for  opaque  skin 


hTb 


= 0 for  diathermanous  skin 


(A  19) 


= (T^  - Tb)/Tb 

0O  j = fictitious  dimensionless  temperature 

(iii)  ®Ni+N2+N3+i,j  = 0 j = 2,  3,  4 ...  (A20) 
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(iv) 


®Nl+\,j  ®N1J  _ ®N1  +2J  ®Nl+l,j 


AX 


N1 


AX 


N1+ 1 


(A21) 


®Nl+N2+l,j  ®N1+N2j 
+N2 


h ®N1+N2+2J  ®NJ  +N2+IJ 

k2  AXN1+N2+ 1 


(A22) 


j = 2,  3,  4 ... 

Where  k1?  k2  and  k3  are  the  thermal  conductivities  of  the  epidermis,  dermis  and 
subcutaneous  tissue,  respectively. 

Since  the  temperature  at  the  (Nl+N2+N3+l)-th  point  is  given  in  Eq.  (A20)  and  the 
fictitious  temperature  at  the  0-th  point  is  given  in  Eq.  (A  19),  rearrangements  of  the  finite 
difference  form  of  Eq.  (A17)  for  i = 1,2,  ...,  Nl,  Nl+2,  ...,  N1+N2,  N1+N2+2,  N1+N2+N3; 

Eq.  (A21)  for  i = Nl+1  and  Eq.  (A22)  for  i = N1+N2+1  provide  (N1+N2+N3)  coupled  algebraic 
equations  for  the  (N1+N2+N3)  unknown  nodal-point  temperature  at  time  j+1.  The  matrix  form 
of  these  equations  is  given  in  Eq.  (A23).  The  elements  of  the  tri-diagonal  matrix  of  Eq.  (A23) 
are  defined  in  Eqs.  (A23a),  (A23b)  and  (A23c).  While  the  elements  in  the  vector  matrix  on  the 
right  hand  side  of  Eq.  (A23)  are  given  in  Eq.  (A23d). 


for  i = 1 
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(2 -bfaj  + fljOjj  + 2710^  + />///(/)  + -i^idJ) 

for  i = 1 

-fl/0/-lJ  + - C/0/+lJ  + F1 VJ) 

for  i = 2,3 TV/ 

0 /or  / = TV/  + 1 


1J 


0 


[4-^]0/j  - c/0/+1j  + F2GJ) 

for  i = TV/  +2,  TV/  +3,...,7V7  +TV2 
for  i = TV/+TV2+1 


(A23d) 


"a/e/-ij  + [4“M0ij  ~ c/e/+w  + ^3ft/) 

for  i = A/  +A2  +2, A/  +A2  +3,..., A/  +A2  +A3-1 

“aN7-^2+N30N7+N2+N3-l,;  +t4  “fyw-t7V2+Atfl  0N7  +N2+N3J  ~ CN1  *m  +N3^N1  +N2  +N3  +1 J 
+ 0A/7  +N2  +N3  +!,/+! ) +F3(N1 +N2+N3,f)  for  i = N1+N2+N3 
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where 


a*  = 


k ' 


vCP, 


At  = 


G N 


Cp 


(Cp)6Ar 


Jk 


hAt 

(Cp)iAxi 


PHI(j)  = 


At 

{Cp)lAx{Tb 

0 


(4;+<W 


for  opaque  skin 
for  diathermanous  skin 


Fk(fJ)  = 


A?/; 

(cp)*r4 


ta/+<W 


0 


for  diathermanous  skin 
for  opaque  skin 


and  fx  is  given  in  Eq.  (All)  which  is: 


Y, 

for  1 < / < N1  + 1 

y2 

for  i > N1+ 1 

The  subscript  k for  a,  a and  F denotes  the  skin  layer  in  the  following  manner: 


for  epidermis 
for  dermis 

for  subcutaneous  tissue 


and  the  subscript  b denotes  the  blood.  Also,  Ax  and  At  in  the  above  equations  are  no  longer  the 
dimensionless  AX  and  Ax,  although  Eq.  (A23)  is  a set  of  non-dimensional  equations. 


62 


Equation  (A23)  can  be  solved  for  the  dimensionless  temperatures  at  time  j+1  using  the 
Thomas- Algorithm,  provided  that  the  temperatures  at  time  j are  known.  In  this  model,  the  initial 
temperature  gradient  is  assumed  to  be  linear  from  the  skin  surface  (32.5°C)  to  the  deepest  portion 
of  the  subcutaneous  region  (37  °C).  Then,  the  new  temperatures  at  j = 2 can  be  calculated. 
Repeatedly,  the  temperatures  for  each  time-step  following  are  found.  The  time-step  employed 
in  this  model  is  1 second.  The  numbers  of  grids,  the  stretching  ratios,  the  thermal  properties,  the 
blood  perfusion  rates  and  the  thicknesses  for  the  epidermal,  dermal  and  subcutaneous  layers  used 
in  this  model  are  summarized  in  Table  A2.  The  surface  heat  loss  coefficient,  h,  given  in  [A3] 
is  4.5  x 10"4  cal/cm2-  sec-°C.  If  the  skin  surface  is  covered  with  garments,  it  is  suggested  to 
ignore  the  heat  loss  coefficient.  Furthermore,  the  emissivity  of  human  skin  is  assumed  to  be  one. 

TABLE  A2  Parameters  Used  in  the  Model 


Parameters 

Epidermis 

Dermis 

Subcutaneous 

Tissue 

Blood 

Number  of  Grid 

6 

24 

15 

- 

Stretching  Ratio 

1.25 

1.25 

1.25 

- 

Overall  Thickness  of  Layer  (cm) 

0.008 

0.2 

1.0 

- 

Conductivity,  k 
(cal/cm -°C*  sec) 

6.1  x {O'4 

12.5  x 10-4 

4.0  x 10"4 

- 

Volumetric  cap. 
Cp(cal/cm3-°C) 

1.032 

0.924 

0.660 

0.954 

Blood  perfusion  rate,  G 
(cm3/sec  per  cm3  of  tissue) 

0.0 

0.00125 

0.00125 

- 

Absorption  coefficient,  y (cm-1) 

180.0 

7.5 

7.5 

- 

In  order  to  test  the  computational  accuracy  of  the  model  and  the  effects  of  the  numbers 
of  grids,  the  stretching  ratios  and  the  time-step,  a simplified  problem  was  studied.  This  problem 
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considers  an  opaque  skin  with  a semi-infinite  boundary.  A constant  heat  flux  of  0.282  cal/sec- 
cm2  and  a uniform  initial  temperature  of  32.5°C  were  used.  The  thermal  properties  of  the  skin 
were  simplified  to: 

kj  = k2  = k3  = 8.83  x 10"4  cal/cm-00  sec 
(Cp)j  = (Cp)2  = (Cp)3  = 0.924  cal/cm3-°C 
Gj  = G2  = G3  = 0 

and  the  surface  heat  loss  coefficient  was  set  to  zero.  Thus  an  analytical  solution  for  this  simple 
problem  can  be  obtained  from  Eq.  (A7).  The  analytical  solution  was  then  compared  with  the 
finite  difference  solutions  for  various  numbers  of  grids,  stretching  ratios  and  time-steps.  The 
different  sets  of  parameters  used  are  defined  in  Table  A3,  and  the  results  of  the  temperature  at 

0. 008  cm  for  the  different  cases  are  shown  in  Table  A4. 

Comparison  of  Cases  a,  e and  k in  Table  A4  indicates  that  for  the  range  of  grid  sizes 
studied,  increase  in  number  of  grids  does  not  obviously  improve  the  accuracy.  However,  change 
of  the  stretching  ratios  does  change  the  finite  difference  solution.  As  indicated  in  Cases  e to  i, 
the  finite  difference  solution  slightly  approaches  the  analytical  solution  when  the  stretching  ratios 
are  increased  until  they  reach  a value  of  about  1.3,  above  which  accuracy  decreases.  A value  of 
1.25  for  the  stretching  ratios  appears  to  work  satisfactorily  for  almost  all  cases  (Cases  b,  c,  g,  j, 

1,  m and  n).  Also,  Cases  b,  c,  d,  e,  1 and  m indicate  that  decreasing  the  time-step  by  half  does 
not  show  obvious  improvement  on  the  accuracy  for  the  heat  flux  of  0.282  cal/cm2*  sec. 

According  to  the  above  analysis,  the  values  of  the  numbers  of  grids  and  the  stretching 
ratios  shown  in  Table  A2,  and  a time-step  of  1 sec  appear  to  be  appropriate  for  solving  this  type 
of  transient  temperature  problems.  The  finite  difference  approach  developed  in  this  section  is 
useful  for  solving  the  skin  model,  and  the  computational  accuracy  is  reasonably  good. 
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TABLE  A3  The  numbers  of  grids,  stretching  ratios  and  time- steps 
used  for  testing  the  model 


Case 

N1 

N2 

N3 

R1 

R2 

R3 

At  (sec) 

a 

50 

50 

20 

1 

1 

1 

1 

b 

50 

50 

20 

1.25 

1.25 

1.25 

1 

c 

50 

50 

20 

1.25 

1.25 

1.25 

0.5 

d 

12 

36 

24 

1 

1 

1 

0.5 

e 

12 

36 

24 

1 

1 

1 

1 

f 

12 

36 

24 

1.1 

1.1 

1.1 

1 

g 

12 

36 

24 

L25 

1.25 

1.25 

1 

h 

12 

36 

24 

1.5 

1.5 

1.5 

1 

i 

12 

36 

24 

1.65 

1.65 

1.65 

1 

j 

8 

28 

18 

1.25 

1.25 

1.25 

1 

k 

6 

24 

15 

1 

1 

1 

1 

1 

6 

24 

15 

1.25 

1.25 

1.25 

1 

m 

6 

24 

15 

1.25 

1.25 

1.25 

0.5 

n 

2 

16 

9 

1.25 

1.25 

1.25 

1 

TABLE  A4  Comparison  of  the  Finite  Difference  and  Analytical  Solutions  for  A Transient  Temperature  Problem 
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BURN  MODEL 

Bum  damage  is  assumed  to  follow  a first  order  reaction  rate  process  [Al],  the  Arrhenius 
equation: 

/ -AE 

Q.  = p J e RT  dt 


where: 

Q 


P 

AE/R 

T 

t 


a quantitative  measure  of  bum  damage  at  the  basal  layer  or  at  any  depths  in  the 
dermis 


= frequency  factor,  sec'1 

= activation  energy  (AE  in  cal/mole)  divided  by  the  gas  constant  (R  in  cal/mole-  K) 
= absolute  temperature  at  the  basal  layer  or  at  any  depths  in  the  dermis,  K 
= total  time  for  which  T is  above  317.15  K (44°C),  the  minimum  injurious  temp. 
Several  bum  criteria  from  the  literature  are  available  for  use  with  this  model.  The  differences 
in  the  criteria  are  mainly  the  different  values  of  P and  AE/R.  These  criteria  include: 

1.  Henrique’s  Burn  Criterion  [Al] 

Based  on  experiments  of  epidermal  bums  and  extended  to  predict  dermal  bum  damage. 
Considers  only  the  damage  incurred  during  the  heating  up  period  of  a bum. 


P = 3.1  x 1098  sec'1 
AE/R  = 75,000.0  K 

2.  Stoll’s  Burn  Criteria  [A2] 

Based  on  experimental  epidermal  bum  data  and  extended  to  predict  dermal  bum 
damage.  It  considers  the  total  damage  incurred  during  both  the  heating  and  cooling  periods. 
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P = 2.185  x 10124  sec-1 
A E/R  = 93,534.9 


for  T < 50°C 


P = 1.823  x 1051  sec"1 
A£//?  = 39,109.8  £ 


/or  7 > 50°C 


3.  M.I.T.’s  Bum  Criteria  [A3] 

Based  on  both  epidermal  and  dermal  experimental  bum  data.  It  accounts  for  damage 
during  both  heating  and  cooling  periods  and  includes  body  fluid  convection  effects. 


P = 1.43  x 1072  sec  1 ► 
A E/R  = 55,387.71  K 


for  epidermis 


4. 


P = 2.86  x 1069  sec'1  - 
A E/R  = 55,387.71  K j 


for  dermis 


Takata’s  Criteria  [A4] 

These  criteria  were  developed  for  dermal  burn  damage  assessment  only. 


P = 4.32  x 1064  sec"1  > 
A E/R  = 50,000  K 

P = 9.39  x 10104  sec-1  ► 
AE/R  = 80,000  K 


for  T < 50°C 


for  T > 50°C 


Morse,  H.  et  al.  [A5]  recommended  that  the  Stoll  criteria  for  epidermal  bum  prediction 
and  the  Takata  criteria  for  dermal  damage  be  combined  for  the  best  prediction  for  cutaneous 
burns.  This  study  will  use  this  recommendation. 
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In  evaluation  the  models  for  this  study,  the  trapezoidal  rule  is  employed  to  integrate  the 
Arrhenius  equation  numerically  to  obtain  the  damage  function  £2.  Once  £2  is  found,  the  degree 
of  bum  is  determined  according  to  the  following  table. 

Assessment  of  Degree  Bum 

£2  < 0.5  at  0.008  cm  No  bum  damage 

0.5  < £2  < 1.0  at  0.008  cm  1st  degree  bum 

£2  > 1.0  at  0.008  cm  and  £2  < 1.0  at  0.2  2nd  degree  burn 

cm 

£2  > 1.0  at  0.2  cm  3rd  degree  bum 

and  the  depth  of  irreversible  damage  is  defined  as  the  depth  at  which  £2=1. 

Since  the  combined  criteria  of  Stoll  and  Takata  have  been  proven  to  give  the  best 
prediction  of  cutaneous  bums,  these  criteria  were  used  in  the  following  analysis  to  test  the 
combined  accuracy  of  the  skin  and  bum  models.  Different  constant  heat  fluxes  were  input  to  the 
program.  The  time  to  2nd  degree  bum  (i.e.  the  time  at  which  £2  at  0.008  cm  reaches  one)  was 
then  predicted  by  the  program  for  each  of  the  input  heat  flux.  Three  different  versions  of  the 
program  were  studied  and  the  results  were  compared  with  the  measurements  from  Stoll  and 
Quanta  [A9].  The  differences  between  the  three  versions  are  shown  in  Table  A5  and  the  results 
are  plotted  on  log-log  coordinates  in  Fig.  A6. 

In  Table  A5  the  conductivities  used  in  the  old  version  and  New  version  1 were  obtained 
from  Henriques  and  Moritz  [A  10].  The  conductivities  used  in  New  version  2 are  the  averages 
of  those  within  the  three  individual  layers  presented  in  Elkins  and  Thompson  [A8].  As  seen  in 
Fig.  A6,  New  version  2 predicts  the  best  results  of  the  occurrence  of  2nd  degree  bum  (or 
blistering). 
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TABLE  A5  The  Parameters  Used  in  the  Three  Versions 


ki 

k2 

(cal/cm  sec  ° 

k3 

C) 

N1 

N2 

N3 

R1 

R2 

R3 

Old  version 

5 x Iff4 

8.83  x Iff4 

3.83  x 104 

12 

36 

24 

1 

1 

1 

New 
version  1 

5 x 104 

8.83  x 104 

3.83  x Iff4 

6 

24 

15 

1.25 

1.25 

1.25 

New 
version  2 

6.1  x Iff4 

12.5  x Iff4 

4.0  x Iff4 

6 

24 

15 

1.25 

1.25 

1.25 

Since  there  is  no  similar  information  available  from  the  literature  about  the  3rd  degree 
burn,  prediction  of  3rd  degree  burn  cannot  be  justified. 


Tolerance  Time  (Sec) 
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Absorbed  Energy  Rate  (cal/cm2  sec) 


Fig.  A6  Measured  and  predicted  human  skin  tolerance  time  to  absorbed  energy  rate 
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